Extracting the Pair Distribution Function of Liquids and Liquid- Vapor Surfaces by 
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We show that the structure factor S(q) of water can be obtained from X-ray synchrotron exper- 
iments at grazing angle of incidence (in reflection mode) by using a liquid surface diffractometer. 
The corrections used to obtain S(q) self-consistently are described. Applying these corrections to 
scans at different incident beam angles (above the critical angle) collapses the measured intensities 
into a single master curve, without fitting parameters, which within a scale factor yields S(q). Per- 
forming the measurements below the critical angle for total reflectivity yields the structure factor 
of the top most layers of the water /vapor interface. Our results indicate water restructuring at 
the vapor/water interface. We also introduce a new approach to extract g(r), the pair distribution 
function (PDF), by expressing the PDF as a linear sum of Error functions whose parameters are 
refined by applying a non-linear least square fit method. This approach enables a straightforward 
determination of the inherent uncertainties in the PDF. Implications of our results to previously 
measured and theoretical predictions of the PDF are also discussed. 

PACS numbers: 61.10.-p, 61.10.Nz 
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I. INTRODUCTION 

The pursuit after the structure and properties of wa- 
ter is obviously driven by its importance as the medium 
that supports life on Earth. Structural studies were 
initiated from the outset, after the discovery of X-rays 
and neutrons 1 , and have continued to these days, as 
technological advances, such as improved X-ray or neu- 
tron sources and computational capabilities have en- 
abled more precise and refined insights^. High en- 
ergy X-ray diffraction techniques for studying liquids and 
glasses^ have also advanced the field 4 . The major- 
ity of previous and current X-ray studies of bulk wa- 
ter have been performed in transmission mode, where 
an X-ray beam travels through the windows of a con- 
tainer before and after scattering from the sample or 
directly scattering through a flowing cylindrical jet of 
the liquid^^ ' 10 ! 11 ! 12 ! 13 ! 14 ! 15 , Narten and Levy and co- 
workers introduced a Bragg-Brentano type diffractome- 
ter to collect the diffraction pattern from a horizontal 
liquid surface in reflection mode that eliminates sample 
holder absorption and scattering 16 . Herein, we extend 
on the reflection mode by using the grazing angle X- 
ray diffraction (GIXD) mode. To conduct experiments 
in GIXD mode, we take advantage of the liquid surface 
(horizontal) diffractometer, that was first introduced by 
Als-Nielsen and Pershanil. 

The quantity of interest in structural studies is the 
liquid structure factor (or the absolute scattering cross- 
section 18 ) S(q), from which the pair distribution function 
(PDF) g(r) is extracted. Experimentally, however, the 
structure factor is only obtained after accounting for mul- 
tiple corrections and assumptions, which may introduce 
systematic errors into the final results^^. The possibility 
of providing a significant simplification in the measure- 
ment of liquid structure- factors was one of the motiva- 
tions to pursue the structural studies described in the 



present study. Compared with transmission mode, our 
experiments in reflection mode do not require corrections 
such as the ones related to the presence of the container 
or 1/r 2 geometric considerations and corrections due to 
absorption and background scattering become substan- 
tially simpler—. In this study, we determine the struc- 
ture factor of water and compare our results with previ- 
ous experiments. More importantly, as the real part of 
the index of refraction for X-rays is smaller than unity, 
the x-ray beam undergoes total external reflection below 
a critical angle a c , with a finite penetration depth. Thus, 
scans above the critical angle predominantly provide the 
bulk water structure factor, while scans below the crit- 
ical angle probe the structure of the top most layers of 
the liquid surface, providing direct information on the 
restructuring at the water/vapor interface. 



II. DETERMINATION OF S(q) IN REFLECTION 
MODE 



A. Reflection mode setup 

Figure [T] shows the setup used to measure the struc- 
ture factor in reflection GIXD mode from a liquid. The 
incoming beam (propagating along the X-axis) hits the 
flat liquid surface at an angle of incidence a with respect 
to the surface, and the scattered beam is detected at an 
angle (3 with respect to the surface and an azimuthal an- 
gle 20 measured from the X-axis. The scattering vector 
is given by 

q = fcrj (cos j3 cos 28 — cos a, — cos j3 sin 28, sin a + sin (3) . 

(1) 
where fco = 27r/A, and A is the X-ray wave-length. Be- 
cause of the isotropic nature of the system, scattering 



scans are presented as a function of the modulus q = |q| 



q = ko \/2 + 2 sin a sin j3 — 2 cos a cos /3 cos 20, (2) 

which for a — /? = gives the known expression q = 
2fcg sin 9. 




FIG. 1: Setup for measuring bulk structure factor of a liquid 
in reflection mode. The incident beam with a wave-vector 
k; hits the surface at a fixed angle a with respect to the 
liquid surface. For bulk measurements a is kept above the 
critical-angle for total external reflection. The scattered beam 
is collected at an angle j3 with respect to the surface and at 
an angle 29 with respect to the X-axis in the XY-planc. 

We use standard calculations for the liquid structure 
factor 20 . In reflection mode, the absolute scattering cross 
section includes a contribution due to scattering from 
capillary waves^i, however this scattering is negligibly 
small compared to bulk scattering above the critical an- 
gle, and can be neglected. The relation between the 
structure factor and the PDF for bulk water (above the 
critical angle) is 



S(q) = (F 2 ) + {F) 2 J Anpr 2 [g(r) 
b 



l] b ^±dr 



qr 



(3) 



p is the number density of the liquid. (F(q)) is the av- 
erage form factor of an H2O molecule (see Fig. [2]). The 
(F(q) 2 ) is the sum of the coherent and incoherent scat- 
tering from each molecule, 



(F 2 ) = (F) 2 + I CS (q) 



(4) 



where Ics 1S the inelastic Compton scattering. For the 
average form factor and the incoherent Compton scatter- 
ing we use the calculated result given in Ref. [23 which 
does not significantly differ from that given in Ref. [2J. 

The raw measured intensity I(q) depends on variables 
such as the incident angle of the X-ray beam and on the 
X-ray wavelength. Its relation to the structure factor- 
s'^) is given by 



S(q) = C 



I{q)-I B {q) 



I V ef f(q;2e,a,l3)P(q;29,a,P) 



(5) 



where C is a scale factor (determined by our analysis de- 
scribed below) Ig is the background intensity measured 
by lowering the surface below the incident beam and oth- 
erwise conducting the same scan as with the sample. Iq 
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FIG. 2: The coherent form factor (F(q)} 2 (dashed line) and 
the total incoherent Compton scattering (dotted line) as given 
in Ref. 1 23| . The solid line shows the sum of the incoherent and 
coherent scattering (F(q) 2 ). 



is the intensity of the incident beam on the sample which 
is monitored immediately before the sample to account 
for fluctuations in the source, and due to the configura- 
tion of the instrument. V e ff(q;29,a,/3) is the effective 
volume of scattering, and P(q; 29, a, (3) is the polariza- 
tion factor. In reflection mode, the angles a and f3 are 
generally very small (< 3 degrees). Therefore, to a good 
approximation, the polarization factor is given by 

P(q; 29, a, (3) w (1 + cos 2 29)/ 2 - unpolarized 

~ cos 2 29 - polarized, (6) 

The evaluation of the effective volume of scattering is 
described in detail in Appendix A. 



III. A NEW METHOD OF EXTRACTING g(r) 
FROM S(q) FOR BULK WATER 

The aforementioned corrections, to within a scale fac- 
tor, yield S(q) from which the PDF has to be extracted 
by solving Eq. ([3]). Traditionally, the method to ac- 
complish this is by inverse integration of Eq. ([3]) , which 
requires extrapolation to q— values where S(q) has not 
been measured. This extrapolation is not free from tech- 
nical problems, as described, for example, in Ref. [9| Eq. 
(9) and (10). 

Herein, we introduce a new procedure for obtaining 
the PDF (for bulk water) that overcomes the problems 
inherent to the inverse integration. The method works 
as follows; we construct a model function g(r) that is 
generated by parameterized Error functions 



JV+l 



(r) = 5^(G<+i-G i )eri 



2 G w+2- ( 7 ) 



The conditions, g(r) = at the origin and g(r) = 1 for 
large r, imply G\ — and Gn+2 = 1, respectively. The 



method starts with one Error function (N = 0) the pa- 
rameters of which arc refined by a non-linear square fit 
(NLSF) method that minimizes a global quality factor x 2 
to obtain the best fit to the numerically calculated S(q) 
Eq. ([3]). An Error function is added iteratively one by 
one, incrementally increasing N + 1 (and the number of 
parameters). Each additional Error function adds three 
more parameters namely, Gi, r* and <Xj. The number of 
Error functions N used at the end of the process is the 
minimum number necessary to fit the data such that the 
addition of another Error function (with its correspond- 
ing parameters) does not improve the quality factor of 
the fit x 2 - The scale factor C in Eq. (J5J) (dependent on 
incident beam intensity) is also a free parameter that is 
refined in this process, it is dominated by the number 
density of water molecules and the molecular form factor 
(elastic and inelastic) . Alternative techniques to directly 
calculate the PDF from the experimentally determined 
structure factor have been introduced in the past 26 . For 
instance, in the empirical potential structure refinement 
technique (EPSR), the parameters of interatomic poten- 
tial energy function are refined to produce the best fit be- 
tween the simulated and measured structure facto r 26 ! 27 . 
Our approach differs in that it does not rely on any the- 
oretical assumptions. 

This process allows to determine uncertainties of the 
free parameters Gi, r% and <Xj which in turn yield the 
uncertainties in g(r). The spread in the values of Gi, 
Ti and Ui reflects the uncertainties associated in extract- 
ing g(r) from a the error-bars of each point in S(q) and 
the finite g-range in S(q). In other words, the analysis 
of the experimental result does not yield a single PDF, 
but a spread of functions that decode quantitatively the 
inherent uncertainties of the experimental results. 



IV. EXPERIMENTAL DETAILS 

The X-ray scattering experiments were conducted 
on the Ames Laboratory Liquid Surface Diffractome- 
ter at the 6ID-B beamline at the Advanced Photon 
Source at Argonne National Laborator y 28 ! 29 . The highly 
monochromatic beam (16.2 keV with energy resolution, 
AE ~ 2eV), selected by an initial Si double crystal 
monochromator, is deflected onto the liquid surface at 
a specified angle of incidence by a secondary monochro- 
mator (Ge(220) single crystal), which is placed on the 
diffractometer— . The synchrotron X-ray beam is highly 
polarized (« 98%) with the electric field parallel to the 
liquid surface, therefore, the polarization correction (Eq. 
(|6|)) is practically given by P(q; 26, a, 0) — cos 2 (26). 

Ultrapure water (NANOpure, Barnstead; resistivity, 
18.1 Mflcm) was used in the present study. The water 
was contained in a Teflon trough, and a glass plate (area 
10 x 5 cm 2 ) was placed in the trough to form a thin water 
film (« 0.3 mm thick) to reduce the effect of mechani- 
cal agitations on the surface smoothness. To test the 
quality of the surface we routinely checked that the re- 



flection from the surface below the critical angle is nearly 
100%. The trough was encapsulated in an air-tight ther- 
mostated aluminum enclosure (T = 294 K), which was 
continuously purged with a flow of helium gas (bubbled 
through water) during the course of the experiment to 
lower background scattering from air. The height of the 
water surface with respect to the beam center was de- 
termined by a reflection from the surface to better than 
±0.005 mm. 



V. RESULTS AND ANALYSIS 
A. Bulk Water 

Figure GDJA) shows the measured raw intensity versus 
q data at several incident beam angles a after normaliza- 
tion to the incident-beam monitor and after background 
subtraction. It is evident that the raw data is depen- 
dent on the incident beam angle a. Normalization of the 
raw data by the effective volume of scattering V e // Eq. 
(|A2[) and by the polarization collapses all data sets at 
different angles a into a single master curve as shown in 
Figure [3jB). This collapse of the curves does not involve 
any adjustable parameters which provides a very strin- 
gent self-consistent test on the robustness of S(q) (within 
a scale factor). 

To determine the PDF, we apply the NLSF method de- 
scribed in Section III to the master-curve shown in Fig. 
[3l and in the process we also determine the scale factor C 
(Eq. ©) which yields the S(q) shown in Fig.g^A). The 
scale factor determines, without any assumptions, the 
maximum value of S(q) at q w 2.0 A -1 is 63 ± 1 which 
compares well with previously reported values^ ' 14 ! 15 . The 
solid line in Fig. H^A) is obtained from the best-fit pa- 
rameters listed in Tabic U with their uncertainties. The 
dashed line is obtained by small variation of parameters, 
within the determined errors. In general we find that the 
most prominent uncertainty in the determination of the 
PDFs is confined to the region of the first peak (« 2.8 A), 
which is mainly due to the finite q— range of the measure- 
ment. Whereas the goodness of the fit is less sensitive to 
small variations in the height and the width of the first 
peak individually, their correlated value, associated with 
the number of nearest neighbors (NN) to a given water 
molecule is. This number of NN is expressed by the sum 
rule relation 



NN = Airp 



g(r)r dr, 



(8) 



where r m ; n is the location of the first minimum of the 
PDF. Figure 0JB) shows the PDF producing the best fit 
to the data (solid line) and another fit within the allowed 
Errors (dashed lines) given in TablelU Despite the spread 
in the height and width of the first peak, all PDF gave 
within error the same number of nearest neighbors 4.7 
±0.1. Figure0jC) shows calculated S(q) to large q- values 
(beyond our measurement range) using the two model 
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FIG. 3: (A) Scattered intensities versus momentum transfer 
q for various incident-beam angles a at T = 294 K. The back- 
ground at all q values is at the 10 -4 level, in the same units 
as shown in the figure. (B) Same data after normalization by 
the effective volume of scattering V e //. All the data at differ- 
ent incident angles collapse to a single master-curve without 
any fitting parameters. The data are also corrected by the 
polarization factor. This is S(q) up to a scale factor. 



PDFs shown in Fig.[4jB). It should be noted that despite 
slight disagreements about the shape and position of the 
first peak, the present and the above mentioned studies 
all agree about the PDF at distances larger than r m , n . 



TABLE I: Parameters that generate the best-fit calculated 
structure factor using Eq. (0) 
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B. Restructured Water Surface 

Performing the GIXD scans below the critical angle 
provides a pattern that is highly surface sensitive due 
to the finite penetration depth of the evanescence wave 
(the penetration depth at a = 0.064° is w 80 A.), and 
the enhancement by multiple scattering, as predicted by 
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FIG. 4: (A) S(q) (at T = 294 K) obtained after scaling the 
data shown in Fig. [3] (circles), the best fit (solid line) and a 
second fit (dashed line) with different PDF but within the 
uncertainty range. (B) Two PDF's used to calculate the best 
fit the S(q) shown in (A). (C) S(q) calculations extended to 
large q values using the two PDFs shown in (B) showing the 
high q values of S(q) are almost identical. 



the distorted wave Born approximation (DWBA) 28 i 30 ' 31 . 
Figure [5] shows two diffraction patterns above and be- 
low the critical angle for total reflection. The scattering 
from the surface differs from that of the bulk ((a > a c )) 
in several respects. First, as q — ■> the intensity di- 
verges due to diffuse scattering from surface capillary 
waves^i. Second, the main peak of bulk water struc- 
ture factor, at q max ~ 2.0 A -1 , is slightly shifted to 
smaller q values suggestive of a larger intermolecular dis- 
tances at the interface compared to bulk water. This is 
in agreement with recent extended x-ray absorption fine 
structure spectroscopy (EXAFS) measurements of water 



microjets, that show the intermolecular 00 distance is 
5.9% larger than that of bulk water—. Third, the shoul- 
der (q ~ 2.7 A -1 ) is less pronounced compared to that 
of bulk water. 

In the following we attempt to modify the method de- 
scribed in Section Hm for a half- filled space with bulk wa- 
ter to examine whether geometrical effects (truncation 
of electron density at z = 0) can give rise to the ob- 
served differences between the scattering from bulk and 
surface. In the absence of a depth-dependent PDF near 
the surface, we assume that the bulk g(r), obtained in 
the present study, is valid everywhere including the sur- 
face. We argue that the scattering consists of two parts, 
as follows, 



I{a, p, 29) = Icw(a, /?, 20) + I b (a, (3, 20), (9) 



where lew and lb represent capillary-wave diffuse scat- 
tering and bulk scattering contributions, respectively. 
The two terms are derived in detail in Appendix [B] 
Numerically evaluating Eq. ^ and the complementary 
equations in Appendix [Bl with a and b as the sole param- 
eters, and g(r) as defined in Table U (as shown in Fig. 
0JB), we obtain a poor fit to the data (dashed-line in 
Fig. [5]). This shows that, Equation © although ade- 
quately describes the diffuse scattering it does not predict 
the observed shift of peak position at (q ps 2.0 A^ 1 ) or 
the observed change of the feature at q « 2.7 A -1 . This 
implies that the PDF of the top most layers is not the 
same as that of the bulk PDF, and a more refined cal- 
culation considering an anisotropic g(r\\,z) needs to be 
used. 




q (A 1 ) 



FIG. 5: Raw GIXD data after background subtraction above 
and below the critical incident angle for total reflectivity as 
indicated. Solid lines are the best fit as discussed in the text. 
Vertical dashed-dotted lines indicate main peak positions of 
the bulk and surface structure factor. 



VI. DISCUSSION AND SUMMARY 



Comparison with previous determinations of 
the water PDF 



As discussed in the previous section, our determined 
PDF is not a single function but a spread of functions 
differing mainly in the height and width of the first peak 
but preserving the number of nearest neighbors (NN). In 
Fig. [6] we compare two of our PDF models with previous 
results ^ 14 ] 15 . Within error, both results are consistent 
with our determined PDFs. It may be argued that the 
relatively large dispersion on the peak height that follows 
from our results is a consequence of the relatively smaller 
range of g-values measured. We have therefore applied 
the method to compute the PDF described in Section III 
to the data of Hura et al. Ref. [fj. The results, shown 
in Fig. [TT] of Appendix [Cl display a somewhat reduced 
but still quite significant dispersion. In particular, the 
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FIG. 6: A comparison of the PDF's from our study and pre- 
vious studies as indicated. The two PDFs of our study are 
the same as those shown in Fig. |4^B). 

number of NN 4.7 ±0.1 is the same as that obtained by 
Hura et al. (NN = 4.7) but slightly different than that 
obtained by Narten and Levy (NN = 4.4) 

Based on our method for extracting PDFs, we argue 
that it will be a daunting task to reduce the uncertainty 
in the first peak of the PDF as very precise measurements 
to large q— values, with accurate control over systematic 
and other errors is hard to achieve as the expected signal 
is extremely low. Figure HJC) shows the minute differ- 
ences between two possible S(qYs at large q- values, which 



give rise to relatively large differences in the first peak of 
the PDF as shown in Fig.gpB). 



ular dynamics ab-initio calculations 4 ^ report a slightly 
more expanded molecular area at the air-water interface, 
in qualitative agreement with the present study. 



B. Implications for theoretical models of water 

Accurate theoretical models of water are of fundamen- 
tal importance for many problems in the physical prop- 
erties of water, in particular in relation to biological pro- 
cesses. In Ref. 15J, the most popular water models were 
compared in great detail with the experimental results 
for the PDF to assess their validity. The main conclu- 
sion drawn was that the TIP5P model 33 provided the 
most accurate description of the experimental data. 

There are some general implications that follow from 
our analysis in this regard. First, our analysis shows 
that the height and width of the first peak of the water 
PDF does not provide a very stringent test to validate 
theoretical models, as existing experimental results can- 
not accurately resolve these parameters, a point that was 
also noted by Jorgensen and Tirado- Rives 34 . A more 
stringent test is to compare the PDF beyond the first 
minimum, as all recent experimental determinations, in- 
cluding ours, consistently yield almost indistinguishable 
PDFs. In this respect, the very popular TIP3P model 35 
does not compare favorably with experimental data while 
the SPC/Fj 3 ^, the other widely used model, compares 
well, and is more consistent with experimental results. 
A more constraining test on the validity of theoretical 
models at distances lower than the first minimum in the 
PDF is provided by the number of NN, which our result 
places at 4.7 ± 0.1. Both TIP5 and SPC/E are consistent 
with these numbers but TIP3P gives 5.1, significantly 
larger than the experimentally extracted value. Given 
its simplicity, it is remarkable how well SPC/E model 
matches the measured PDF. We conclude that whereas 
experimentally determined water PDFs provide valuable 
guides to test the validity of water models, there is still 
some uncertainty in the experimental region at distances 
smaller than the shell of NN i.e., the first minimum in 
9(r). 



C. Theoretical predictions of the properties of 
interfacial water 



The description of the structure of the air-water in- 
terface by theoretical models is quite challenging. The 
dipole moment of the water molecule strongly depends 
on the environment. In bulk, the dipole moment of a 
water molecule is about 2.4D, while in the vapor phase 
becomes 1.8D. Rigid models, such as the SPC/E, TIP3P 
or TIP5P have a fixed dipole moment, irrespective of 
whether the water molecule is in bulk at the interface or 
in a gas phase. In fact, the surface tension of the vapor- 
water interface calculated with the most popular rigid 
water model a 38 i 39 show significant disagreement (20% or 
more) when compared with experimental results. Molec- 



D. Summary 

The objectives of the present study were twofold: first, 
to demonstrate that the structure factor of liquids can be 
accurately determined by GIXD in reflection mode, both 
to investigate the structure of the bulk or the interface, 
and second, to introduce a new method to compute the 
PDF from the structure factor in a way that allows to 
assess the intrinsic errors associated with the structure 
factor 5(g). 

One advantage of measurements in reflection mode is 
that the raw intensities are very close to the actual 5(g) 
to within a factor, especially at large incident beam an- 
gles, as argued earlier- 1 ^,. The major correction needed is 
the effective volume of scattering, which under suitable 
choices of slits could even be reduced to a trivial scal- 
ing factor independent of q. This should be contrasted 
with transmission measurements, where the measured in- 
tensity is dramatically different from the structure factor 
due to the effect of the container, and geometry (see Fig- 
ure 4 in Ref. |Tj). A unique advantage of the X-ray in 
reflection- GIXD mode is that it can be applied to de- 
termine the structure at the vapor/liquid interface by 
adjusting the incident beam below the critical angle. In 
this configuration the evanescent wave scatters mainly 
from the topmost layers at the surface, and provide valu- 
able information on the restructuring at the gas/liquid 
interface. 

Our experimental results show that the bulk PDF does 
not describe the surface scattering data correctly, sug- 
gesting a restructuring of water molecules at the va- 
por/water interface. In particular, our results show that 
the water molecules at the top most layers are more ex- 
panded. 

In the present study, we have also introduced a 
method based on a non-linear-least-square refinement 
to obtain the PDF from the experimental 5(g). This 
method avoids the problems associated with computing 
the Fourier transform by extrapolations to g-values that 
are not measured, and provides the uncertainties associ- 
ated with the PDF that are compatible with the mea- 
sured experimental data for 5(g). We hope that our 
findings will initiate future experimental and theoretical 
studies of liquids surfaces in general. 



APPENDIX A: EFFECTIVE SCATTERING 
VOLUME, Veff 

The volume of scattering is defined as that region of 
the illuminated sample whose scattered rays are actually 
detected. This volume depends on the angles of the in- 
cident and scattered beam with respect to the surface, 



the apertures of the detector, and the attenuation length 
through the sample (i.e., X-ray energy). In our setup, the 
cross section of the incident beam is rectangular defined 
by two sets of slits with vertical and horizontal opening 
Wi and di, respectively. Typically, Wi ~ 0.01 — 0.1 mm, 
and di ra 1 — 2 mm. The area Ai of the incident beam 
footprint on the liquid surface varies with the angle of 
incidence a 



Ai(z;a) 



w i d l 



dil. 



(Al) 



as illustrated in Fig. \7\ Although the area A t is preserved 



side view 




FIG. 7: (A) Illustration of a side view of the beam footprint. 
As the beam penetrates the bulk the center of the foot print 
is shifted along X. (B) Top view of the beam footprint (z — 0) 
cross section with the footprint of an out-going beam at angle 
29. (C) Same as (B) but for a footprint of the incident beam 
at finite z value. 

as the beam penetrates the bulk of the liquid, the center 



of the illuminated rectangle shifts away from the sample 
axis of rotation, changing the effective volume of scatter- 
ing. The outgoing beam slits, typically Wf ~ 2 mm and 
df ~ di vertically and horizontally, respectively, form a 
footprint that is longer than that formed by the incident 
beam. The effective volume is an integral of the over- 
lap area between the incoming and outgoing footprint at 
each z value weighted by the attenuation length of the 
incident and scattered beams, as follows 



V t 



eff 



A(z, ai 29)e~ z/L dz 



(A2) 



where L = (D(a)D(f3)/[D(a) + D(/3)} is the effective 
attenuation length the beam, and n ~ 8 ensures the con- 
vergence of the numerically calculated integral. The at- 
tenuation length into the bulk at an angle of incidence a 
is given by 



D(a) = 1/Im(k z ) 



where, 



and 



k z (z) = koy sin a — 25 — i2j, 



(A3) 



(A4) 



S = ±- £ N 3 r,\ 2 f' (A) ; 7 = ^ £ N^f- (A) 

3 3 

(A5) 
where Nj is the number density of atom type j with form 
factor / and absorption factor /■ . 
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FIG. 8: Attenuation length for an external incident beam 
(solid line) and for an internal incident beam (dashed line) 
at the vapor /water interface, calculated by Eq. lA3l for a 16.2 
keV X-ray beam. The two curves converge at angles larger 
than the critical angle for total reflection. (Arrow indicates 
the location of the critical angle) 

Figure [8] shows the attenuation length as a function of 
the incident angle of a 16.2 keV X-ray beam propagat- 
ing from the gas phase onto the liquid surface. Below a 



critical incident angle for total external reflection a c the 
beam at the surface is evanescent; it penetrates to a finite 
depth (D(a) < 80 A) into the bulk and emerges almost 
totally without transmission (a small fraction, less than 
0.1%, depending on the surface- roughness, scatters as 
diffuse scattering). Above the critical angle, the beam 
is mostly transmitted but it is also attenuated by the 
absorption coefficient (7) of the liquid. The beam that 
scatters from the bulk at angle /3 emerges through the 
surface with no total reflection, as implied in Figure [51 
The two curves, D(a) and D((3) converge at large angles. 
The effective scattering volume is calculated from 
Eq. (|A2|) by determining the effective scattering area 
A(z, a, 26) at depth z. This area is polygonal in shape 
(see Fig. [7]) and is determined from the geometrical con- 
straints 



1 z I z , . „. 
- ~ <x<- , A6 

2 tana 2 tana 

which account for the shift of the footprint center as the 
beam penetrates into the bulk in the z direction, 



di di 

2 " y ~ 2 



and 



2sin(26») tan(20) 



< x < 



df 



2sin(2(9) tan(2(9)' 



(A7) 
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FIG. 9: Effective volume of scattering using Eq. (|A2[) at vari- 
ous angles of incident beam a and /3 = 0.27°. The calculation 
is for a 16.2 keV beam, and slit parameters di = 1.8, d/ = 2 
mm, w — 0.05 mm. 

Figure [9] shows the effective scattering volume as a 
function of the angle 26 for different values of incident 
beam angle a. When a is relatively large, the footprint 
becomes smaller than the detector aperture and as a re- 
sult, almost all of the illuminated area is detected, that 
is, the effective volume of scattering is almost insensi- 
tive to the scattering angle 26. As a becomes smaller, 
the overlap of incoming and outgoing footprints is more 
complex, resulting in a stronger variation of V e ff on 26, 
as intuitively clear from Fig. [7J 



APPENDIX B: CALCULATION OF 

SCATTERING FROM HALF-FILLED SPACE 

SURFACE ASSUMING BULK PDF 

According to the diffuse scattering theor y 21 ' 22 , 

lew (a, P, 26) = aPA(0, a, 26) 
(2q xy + Aqxy) 71 - (2q xy - Aq^ 



2 1 ' 



P 



oc — (q xy > Aq xy ; 77 -> 0), (Bl) 

where a is a scale factor, P is the polarization factor, A 
is the effective illumination area at z = (see Fig. 7), 

Qxy = JqJ+~q%, Aq xy = k cos6A6, rj = If^g 2 , and A6 
is the acceptance angle of the detector. 

By assuming g(r) is the same everywhere including 
the region close to the surface, the bulk scattering from 
a half-filled space with a penetration depth L (see Eq. 

SJ) can be written as 

I b (a,0,26)=bP A{z,a,26)w(q,z)dz, (B2) 

Jo 



where & is a scale factor, n = 8 as discussed in Eq. (|A2|) . 
and w(q, z) is the scattering intensity from one water 
molecule at a depth z with respect to the surface given 
by 

/oo 
-z 
xe -(2z+z')/2L dz , 



r a i-2-k 

< 1 dr\\ 
/o jo 



d6pr\\(g{Jrj + z' 2 )-l)c 



-iq z z 



-iq xy r\\ cos 9 



(B3) 

where p is the number density of the water, e~ z ' L and 
e -(2z+z )/2L are ^- ne attenuation factors for the inelastic 
and elastic scattering, respectively. To make the numeri- 
cal integration more efficient, we define a sphere of radius 
a (~ 12 A) outside which g(r) = 1 (Fig. OH for L > a . 
One should note that w(q,z) = e~ z/L S(q) (see Eq. ©) 
for z > do- For large incident angles, L>ao, Eq. (|B2|) 
can be simplified to 

I b ~ bS(q)P / e~ z/L A(z, a, 26)dz __ bS(q)PV e 



eff> 

Ja 

(B4) 
which is practically the same as Eq. ([5]). For small in- 
cident angles, q z ~ so, to simplify the integrals in Eq. 
(|B3|) we assume q z — yielding 



w(q,z) 



z/L < F 2 > + < F > 2 

e -{2z + z')/2L dz , 



ra 

x / dm 
Jo 



|27rpr||(ff(^rf + z' 2 ) - l)Jo(wu), (B5) 
where Jq is the zero order Bessel function. 




FIG. 10: Illustrations for attenuation factors and the integra- 
tion range of Eq. IB3I 



APPENDIX C: APPLICATION OF THE NLSF 
METHOD TO OTHER S{q) MEASUREMENTS 
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FIG. 11: (A) S{q) data from Ref.^. (circles). Solid- and 
dashed- lines are the best fits using the method described 
in Section HIT] B) Two extracted PDFs from (A) that fit the 
data equally well. Despite the higher range of the measured 
q-values, possible uncertainty in height, width, and position 
of the first peak in g(r) is evident. 



FigureQTJA) shows two slightly different fits (solid and 
dashed lines) to the S(q) data in Ref . \14 using the method 
described in Section llTTl to extract the PDF. Although the 
two fits are practically of the same quality they produce 
different shapes g(r)'s in particular near the first peak, 
as shown in Fig. [TlTB). We show this to demonstrate 
that even with data measured to larger q values than in 
the present study, some ambiguity in the evaluation of 
the first molecular shell around a water molecule is still 
present. 
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